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Differences between vector potentials in different gauges contain no dynamics in both classical and 
quantum electrodynamics and chromodynamics. Consequently, once gauge invariance is established, 
results calculated in non-covariant gauges can be expected to agree with results obtained in covariant 
gauges in all Lorentz frames. We show in particular that canonical quantization in the Coulomb 
gauge can be used without giving up explicit gauge invariance. Quantization in the Coulomb gauge 
is particularly simple because it involves only the two transverse photons/gluons present in all 
gauges. These transverse photons/gluons reside on a 2-dimensional physical plane in momentum 
space perpendicular to the photon/gluon momentum k. Explicit expressions are given for the basic 
momentum, spin and orbital angular momentum field operators of photons/gluons in the Coulomb 
gauge. Their properties are discussed in some detail. In particular, these field operators are shown 
to be more complicated than the corresponding operators in quantum mechanics which they also 
contain. 
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I. INTRODUCTION 

In classical electrodynamics (CED), the electromag- 
netic (e.m.) fields are uniquely determined by the 
Maxwell equations. The Maxwell theory contains a hid- 
den Lorentz symmetry first uncovered by Einstein f]] 
which can be displayed explicitly by using a vector poten- 
tial. Unfortunately, the vector potential is not uniquely 
determined by the Maxwell equations. There are in- 
finitely many vector potentials or gauges that generate 
the same e.m. fields [3] . This redundant gauge degree of 
freedom has to be removed by a choice of gauge before 
any actual calculation using the vector potential. It is 
obvious, however, that the e.m. fields generated by any 
two different vector potentials Al^{x),n — 1,2, must be 
the same, so that the difference A'^^{x) = ^2 (x) — A^l{x) 
is a function of the spacetime location x that carries no 
dynamical information whatsoever. Such a nondynami- 
cal (hence the subscript "nd") vector potential has the 
functional form 9^A(x), where A(a;) is any twice dif- 
ferentiable local function of x, because the e.m. field 
F^^ = df^d^Aix) - d^df^Aix) generated by it is identi- 
cally zero. In this paper, the word "nondynamical" will 
be used to describe a quantity that contains no dynamics. 

Quantum ED (or QED) differs from CED by including 
the quantum properties of matter and radiation. The 
charged particle is now an electron described by a rel- 
ativistic Dirac wave function -0 spread out in space and 
accessible through a kind of square-root decomposition of 
the charge density: p{x) — etp^ {x)ip{x), where e = — |e| 
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is the electron's charge. The electron has the quantum 
momentum density distribution in space ip'^p^^, where 
the quantum (or quantum mechanical) momentum oper- 
ator = /i has to be placed between the two wave 
functions. This momentum density does not change if a 
global phase factor exp(iA) with an x-independent phase 
A is added to ip. On the other hand, if the phase is a twice 
differentiable local function of x, the momentum density 
changes to tp'^{p^ + d^X)ip. Thus the expression is not 
invariant under a local phase transformation described 
by a phase function X{x) that contains no dynamics. 

In QED, the vector potential plays a role more funda- 
mental than the e.m. fields because it has the Lorentz 
vector structure of spacetime itself. The Lorentz force 
law of CED shows that the charge-field interaction 
changes the point particle's canonical momentum^'^ into 
the dynamical (usually called kinematical or mv [3] ) mo- 
mentum p^ — eA'^(x). In QED too, dynamics now ex- 
pressed as an electron-photon (ej) interaction can be 
included explicitly by using the dynamical momentum 
density distribution tp^{p^ — eA^)tp in space. Under a 
simultaneous gauge transformation of the vector poten- 
tial and a local phase transformation of the wave func- 
tion, the operator part of the expression is changed to 

_|_ giJ-\ _ e(j4'^ -I- d^A). The original expression is re- 
stored if the two changes are chosen to cancel each other 
perfectly: d'^X = ed'^A 0, This is the archetype 

of manifestly gauge invariant expressions usually used in 
QED, including the momentum and angular momentum 
(AM) densities proposed by Ji 0] and by Chen and Wang 
[sl in the nucleon spin problem. With manifest gauge in- 
variance, one is assured that calculations made in any 
gauge give the same physical results as those obtained in 
any other gauge. 
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However, neither the nondynamical phase function 
X{x) nor the nondynamical gauge function A(a;) contains 
any 67 interaction dynamics. The dynamics is already 
present in the Lorentz force law. We are only concerned 
with its invariance under a nondynamical gauge trans- 
formation, so that this gauge transformation itself does 
not add extra dynamics into the problem. Such a re- 
quirement has the same significance whether or not the 
A that appears in the expression contains any dynamics. 
Moreover, the nondynamical expressions are closer to the 
expressions obtained directly from the QED Lagrangian. 
For these reasons, we have previously proposed leaving 
out the 67 interaction altogether in constructing a gauge 
invariant version Q of the operators proposed by JafFe 
and Manohar , where the gluon orbital AM and spin 
operators are separated. In contrast, Wakamatsu |ll| has 
more recently advocated the dynamical version of these 
operators. 

A somewhat more general conceptual scheme for real- 
izing nondynamical gauge invariance works as follows: 
Let a general vector potential be written as = 



quan 



^nd' where A^^^^^ is the vector potential in the 



gauge where second quantization will be realized. Then 
the nondynamical momentum density ij}''{p'^ — eA'^^)ip is 
manifestly gauge invariant. Moreover, the nondynamical 
quantum momentum operator — eAj^^ satisfies all the 
commutation relations of alone. In contrast, differ- 
ent components of the dynamical momentum — eA^ 
do not commute with one another Hence our non- 
dynamical, gauge-invariant QED momentum operator is 
also a canonical momentum. 

We note in passing that in this paper, we call the quan- 
tum operator p = V/i a "nondynamical" operator in the 
sense that its eigenfunctions exp (jk • x) are momentum 
eigenfunctions describing a free particle of momentum 
k in free space. It is technically a canonical momen- 
tum. Feynman calls p simply and unambiguously the 
"p- momentum" [sj. However, it is sometimes also called 
perhaps confusingly a "dynamical momentum" Q , a us- 
age we do not follow. In this paper, a "dynamical" op- 
erator refers to one containing dynamics or interactions, 
unlike a "nondynamical" operator. 

The nondynamical phase factors considered here are 
functions only of spacetime locations. Their phase func- 
tions are obtained from a phase integral that is path- 
independent or integrable. There is another type of phase 
factors involving path-dependent or nonintegrable phase 
integrals that can add dynamics to the fermion field 
0, [3 • It is amusing to note that in the local gauge 
theory of interactions [61] , local gauge invariance stated in 
terms of the nondynamical, path-independent phase fac- 
tor is used to deduce the mathematical structure of the 
dynamical, path-dependent phase factor. However, the 
local gauge theory of interactions is not an immediate 
concern of this paper. 

The purpose of this paper is to present and discuss the 
properties of the nondynamical, gauge-invariant opera- 
tors in QED and QCD 9]. Our method is very fiexible 



because it does not exclude the use of dynamical oper- 
ators. On the other hand, the separation of the non- 
dynamical momentum from the 67 interaction seems to 
provide a more practical calculation scheme from the per- 
spectives of experimental measurement, AM algebra and 
perturbation theory. Our presentation will be detailed 
and self-contained, in order to appeal to a wide reader- 
ship. A historical note: According to Jackson and Okun 
the Lorenz gauge, usually called the Lorentz gauge 
and attributed to H. A. Lorentz, was first proposed by 
L. V. Lorenz in 1867. 

We shall begin in Sect. II with QED. Its natural habi- 
tat is the four-momentum space where the Maxwell equa- 
tions appear as coupled algebraic equations that can be 
solved analytically. Explicit formula for the vector po- 
tential will be given for a family of gauge choices that 
include the Lorenz and Coulomb gauges. These formulas 
will be needed to construct the more complicated opera- 
tors in QCD. 

None of the vector potentials At^ of the family is a 
Lorentz vector except those in covariant gauges such as 
the Lorenz gauge. (Note that A^ is a Lorentz vector if it 
satisfies the covariant gauge condition C^A^ — Lorentz 
scalar, where C'^ is a known Lorentz vector such as 
or x^.) Hence manifest Lorentz covariance and gauge 
invariance can coexist only for canonical quantization in 
a covariant gauge such as the Lorenz gauge. The trouble 
is that the covariant A'^ then generates four covariant 
photons: Among these photons, the longitudinal photon 
is unphysical and the time-like photon is virtual, never 
free of its source charge. Only the two remaining photons 
are transverse photons that can propagate in free space. 

In all non-covariant gauges, manifest Lorentz covari- 
ance of the vector potential is lost. Photon quantization 
must then proceed sequentially by first using a covari- 
ant gauge to ensure Lorentz covariance of the vector po- 
tential in different Lorentz frames. Only then can one 
work in a single Lorentz frame to construct manifestly 
gauge invariant operators. From this perspective, gauge 
invariance implies Lorentz covariance, but Lorentz co- 
variance does not imply gauge invariance. That is, gauge 
invariance in one frame includes a covariant gauge where 
Lorentz covariance can be restored. On the other hand, 
Lorentz covariance in different frames alone does not in- 
volve gauge invariance at all. 

The 3-vector potential is particularly simple in the 
Coulomb gauge where it contains only the 2-vector part 
Ac ~ 'A± — ^phys in momentum space perpendicu- 
lar to the propagation vector k of the light wave. This 
2-vector potential is common to all the gauges in our 
chosen family of gauges. It contains all the physics of 
the two transverse free photons. Consequently, we shall 
call this transverse plane the "physical plane" (hence the 
subscript "phys" .) Although ^phys is defined only after a 
Lorentz frame has been chosen, it nevertheless develops 
a hidden Lorentz symmetry under the special circum- 
stances of QED, namely that .4^hys is a Lorentz scalar 
for all the ^phys of different Lorentz frames. 
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The construction of momentum and AM operators in 
QED using Noether's theorem directly in terms of the 
vector potential as a field variable gives results where 
the nondynamical electron and photon terms are cleanly 
separated. Dynamical operators can next be introduced 
if desired by adding a canceling action-reaction pair of 67 
interaction terms and associating one partner term with 
the electron operator and its canceling partner with the 
photon operator. 

Our nondynamical photon orbital AM and spin oper- 
ators are well-defined gauge-invariant operators. Each 
operator has a well-defined component along the photon 
momentum k that describes the effect of 2-dimensional 
rotations about k. They are therefore separately useful 
for problems with cylindrical symmetry about k. The 
physical transverse photon states can be expressed in 
terms of eigenstates of the total AM J = L + S. We shall 
explain why the absence of longitudinal photon states in 
free space makes it impossible to observe the photon L 
and S separately and completely in free space. 

The extension of the above treatment to QCD is con- 
ceptually the same, but technical more complicated be- 
cause of the nonlinear interactions among the vector po- 
tentials. We show why the family of gauge conditions 
used for CED can also be defined in CCD. Then an 
explicit perturbative solution in momentum space for 
the nondynamical part A'^^ can be obtained. (Note 
that ^{^ui-e @) Ell is j^st -^nd special case 

where quantization is realized in the Coulomb gauge.) 
In particular, the vector potential Aj^ transverse to the 
gluon propagation direction k is found to be gauge- 
independent, just like CED, at least for our perturba- 
tive solutions. We therefore do not use the conceptual 
method for QCD proposed in [1, [H, [3, but construct 
QCD momentum and AM operators directly in complete 
analogy to the QED operators. Canonical quantization 
on the physical plane in the Coulomb gauge shows that 
the QCD operators, though nominally nondynamical like 
their QED counterparts from the perspective of gauge 
transformations, nevertheless contain nonlinear interac- 
tion terms proportional to the color coupling constant 
g. Hence the complete removal of dynamical effects from 
the QCD operators seems impossible. 

Sect. IV contains concluding remarks, where we touch 
upon the recent criticism of our operators published by 

ji 



obtained in Sect. E. The units h 



1 are used. 



A. Vector potentials in classical ED 

In a chosen Lorentz frame, the vector potential has a 
particularly simple structure in the 4-dimensional Fourier 
space kt" = (w,k) [ll,^: ^^(fc) = -l-^j.). It 

contains a two-dimensional part transverse to the 
propagation vector k and known to be gauge indepen- 
dent because it is completely determined by the gauge- 
independent magnetic induction: 



rk X 



'_L- 



(1) 



Once known, A.± determines £± from the Faraday effect: 



£± = iLuA._ 



(2) 



There is now only one e.m. field 5|| = £||ek left to be 
related to the two leftover components A'^ and A\\ . The 
remaining relation between them is usually written as 



£11 = -i\k\A° + iLoAii, 



(3) 



showing that there are infinitely many choices of A'^ and 
^11 for the same f || . This is the redundant gauge degree 
of freedom. 

We are interested in a family of vector potentials de- 
fined in any Lorentz frame by the following real linear 
relation, to be called the generalized velocity gauge con- 
dition: 



A 



{■"a) 



k 



±1, 



(4) 



where the gauge parameter a — is expressed in 

terms of a gauge velocity Vg that is a Lorentz scalar. The 
choice s = 1 has been studied in (T9l - [2l| . We are here 
extending the family by including the subfamily with the 
signature s = — 1. (If this s = — 1 signature factor is 
absorbed into Vg, the gauge velocity will become purely 
imaginary.) The Coulomb gauge ac = appears in both 
subfamilies. The gauge independence of in Eq. ([3]) 
requires that [l^ 



rA%. 



(5) 



II. GAUGE INVARIANT OPERATORS IN QED 

We begin in Sect. A with classical ED. The properties 
of vector potentials and their explicit functional form in 
a certain family of gauge choices will be derived. In Sect. 
B, we construct gauge invariant nondynamical and dy- 
namical momentum operators and briefly discuss their 
merits. Angular momentum operators are constructed in 
Sect. C. Photon AM operators are examined in detail 
in Sect. D, and their eigenstates of good total AM are 



That is, Eq. ([5]) is a solution of the simultaneous linear 
equations ([3l |4]) for the unknowns ^'^(""s) and ^y"^' in 
any Lorentz frame. 

Among these gauges, the Coulomb gauge provides the 
most economical description because ^c|| = and Aq = 
f||/(ia;). That is, Ac = A± — Aphys contains only the 
gauge-independent and physics-containing part common 
to all gauges that resides on the physical plane. We shall 
therefore call the Coulomb gauge the physical gauge. The 
vector potential in any gauge is then — A'^^^y^ + A!^^^^, 
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"pure" refers to the nondynamical 



where the subscript 
pure-gauge part. 

It is interesting to note here that in QED, A.± on the 
physical plane develops a hidden Lorentz symmetry: 
is a Lorentz scalar for one-photon states. This symmetry 
appears because in QED, A'' in a covariant gauge is to be 
quantized into a one-photon operator. Consistent with 
this objective, its Fourier space (w,k) is specialized into 
the 4-momentum space of a massless photon satisfying 
the condition w = |k|. For example, the Lorenz gauge 
condition then reads 

Am = Al (6) 

Then the Lorentz scalar Aj^ of the vector potential in the 
Lorenz gauge is just 



Al 



Al 



Al 

\ 2 



AlA 



{A 



-Ai 



(7) 



[We use the metric g'^'^ — diag(l, —1, —1, —1) for covari- 
ant quantities. For 3-dimensional vectors, however, we 
do not use a tensor notation so that spatial indices have 
the same meaning as superscripts or subscripts. Their 
location is determined only by the readability of the ex- 
pression.] Under Lorentz transformations, both the ori- 
entation of the physical plane as defined by the propa- 
gation direction et and the actual direction of Jl± on 
the physical plane can change, but not the length A± 
contained in A'^ that is common to all the gauges under 
consideration. 

Furthermore, starting with the Lorentz vector poten- 
tial A'^ in a new Lorentz frame, we can gauge trans- 
form in a standard way to any other velocity gauge [l9| . 
The resulting vector potential in this new gauge in the 
new Lorentz frame must satisfy the gauge condition of 
Eq. dH), by definition. In the Coulomb-gauge, for exam- 
ple, Ac\\ = and therefore the Coulomb-gauge condition 
k - ^ = are guaranteed in the new Lorentz frame. From 
this property in classical ED, one can expect that QED 
formulated in any velocity gauge gives the same physical 
results even though the vector potential is in general not 
a Lorentz vector. For an explicit demonstration of the 
Lorentz invariance of QED in the Coulomb gauge, see 

[Hill. 

In the s — —1 subfamily of velocity gauges, the choice 
a = I gives rise to the covariant light-cone (Ic) gauge 
where Aic\\ = —A'l^ for one-photon states. Consequently, 
Af^ = —A]_ holds also. 

Returning to CED, we note that another hidden 
Lorentz symmetry can be said to exist in the A'^^ be- 
tween any two gauges, namely it generates no e.m. field 
i-^nd ~ ^) Lorentz frame even though it is not a 

Lorentz vector. 

The family of vector potentials even in one gauge on 
one Lorentz frame is very large, for the gauge condition 
(HI can be written more generally as 



A\^^^+X^sa^iA°^^^^ + ^X 



where X: a solution of the "wave" equation 

(|k|2-W)x = 0, (9) 

defines a subfamily of restricted gauge transformations 
in the same gauge. The gauge velocity Vg contained in 
the gauge parameter a = is the nonphysical wave 
velocity appearing in Eq. 

Because the two components ^^("s) and A^^^'^ are de- 
scribed by a nonphysical wave velocity Vg, the resulting 
j[i^{vg) jg j^Q^ Lorentz vector except when Vg = 1, in 
the Lorenz gauge when s = 1 and in the covariant light- 
cone gauge when s — —1. However, The two remaining 
components on the physical plane are gauge indepen- 
dent and satisfy physical wave equations involving c = 1. 
This is the mathematical reason why ^± contains hid- 
den Lorentz symmetry even though it is not a Lorentz 
4- vector. 

The Lorenz and covariant light-cone gauges are also 
members of the generalized covariant gauges satisfying 
the gauge condition C^A'' = Lorentz scalar, where A'^ 
is a Lorentz vector: Here is a known Lorentz vector 
and the word "generalized" refers to the use of a Lorentz 
scalar that is not zero. Most covariant gauges of practi- 
cal importance use a zero Lorentz scalar. These include 
the Lorenz gauge [C^ = = (w,— k)], the covariant 
light-cone (Ic) gauge [C^ = fc^ = (i^jk)], and the rela- 
tivistic Poincare gauge [C^ = — (t, — x)]. Covariant 
gauges like the relativistic Poincare gauge are defined in 
spacetime where the physical plane in the Fourier space 
k cannot be isolated explicitly. That is, only the Fourier 
transform lies on a 2-dimensional plane in the Fourier 
space k. 

The relation between A\\ and A'~' can also be nonlinear: 
^11 — /(A*^), where f{x) is a nonlinear function of x. 
The main complication of such a choice is that the gauge 
independence condition of Eq. ([3|) in general has multiple 
roots. 

Last but not least, the classical e.m. momentum and 
AM densities can be written in the following useful dif- 
ferential triple vector product forms when expressed in 
terms of the vector potential A: 



E X (V X A) = 



X X [E X (V X A)] 



(10) 



E^VA^ + (V • E)A 
E^(-K X V)A3 

+ (V-E)(xx A)+Ex A 
- V^'[£;^'(x X A)]. (11 



(8) 



Equation (fTTj) is easily derived from Eq. (|T0)) . These 
elementary classical expressions already contain the dif- 
ferential operator and gauge-dependent substructures we 
shall find in QED. Most presentations of the subject use 
these identities to change the classical expressions di- 
rectly to their QED relations. We shall work directly 
with the QED Lagrangian, however, in order to provide 
a different perspective. 
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B. Gauge invariant nondynamical and dynamical 
momentum operators in QED 

The QED Lagrangian density is 

Ff,, = df,A,~d,A^, (12) 

where the field variables are tp-,ip and A^. According 
to Noether's theorem, the invariance of C under trans- 
lation generates a conserved energy-momentum tensor 
T^"^. The e and 7 momentum densities 



r 



On 



rpOn I rpOn 



(13) 



contain no direct contribution from the 67 interaction 
term in £. Here n, j are Cartesian coordinate labels in (3- 
dimensional) space. Repeated indices imply summation. 

Since the classical Maxwell theory is Lorentz covari- 
ant, we may next choose a Lorentz frame and work with 
the momentum operator in 3-dimensional space, provided 
that gauge invariance is verified afterwards for all phys- 
ical operators. This is the price paid for bypassing the 
gauge independent B field. The total QED momentum 
operator is then 



= J i^^p^d^x + J[E^{VA^)]no.^d^x, (14) 

where p = V/i is the quantum canonical momentum, 
and 



[E^{VA^)]a 



e„„-i[i?^'(VA^) + (VAJ")i?^] (15) 



is Hermitian. For notational simplicity, the subscript 
"Herm" will not be shown henceforth. The photon mo- 
mentum has a gauge-invariant physical part and a non- 
dynamical, gauge-dependent or pure-gauge part: 

p7 — p7 1 p7 

phys ' pure 

= J E^^{VA:'^)d^x, + J S,|(VA||)d3a;. (16) 

Note that four terms appear on writing E = +E|| and 
A = A_L -l- A||. The two cross terms not shown vanish 
separately: Their structures are more transparent when 
their space integrals are re-written as momentum inte- 
grals. Each integral then involves perpendicular vector 
fields £*i and A.\\ (or ^jj and .4.^) that form a zero scalar 
product when V — )• ik is moved out of their way. 

The gauge-dependent term PJ^m-o can be moved to the 
electron momentum operator to make the latter mani- 
festly gauge invariant. This is done by first express Pp^rc 
in a familiar form with the help of the differential vector 
identity used with E E|| and V x A|| = 0: 



The volume integral over space of the last term on the 
right can be changed into a surface integral "at infinity" 
by using the Gauss theorem (/ d'^xV-' -> Jgda^). The 
surface integral vanishes if the fields are vanishingly small 
on the large surface at infinity. Then 



p7 

pure 



ip^ eAnip d^x 



(18) 



where use has been made of the Gauss law V-E|| = eip'^ip. 
Added to P"^, this nondynamical longitudinal momentum 
gives manifest gauge invariance to the nondynamical elec- 
tron momentum 



pe ^p7 

' pure 



J i^'^ip- eA||)V'd^x. 



(19) 



Once gauge invariance is assured, we can work in the 
Coulomb gauge where the electron momentum is just 
P*^ . Note that the complete gauge invariant 4-momentum 
density is V'^(p^ ~ eA^'m.o)V'i where in momentum space, 
both fc^ and -A^^^^g have only nonzero fi — 0,3 compo- 
nents. Their identical mathematical structures explain 
why they appear naturally together. 

To get Ji's dynamical electron momentum operator: 



dyn 



(20) 



we have to add to our nondynamical PJj^ a transverse 
momentum — P'p' from the ej interaction, where 



P7 = J i/jteA^Vrf^a;. 



(21) 



The canceling partner P'p' of this action-reaction pair of 
terms is then added to the photon momentum Pp^yg to 
give a dynamical photon momentum operator 



p7 _ p7 

dyn phys 



,e7 



(22) 



Such a re-definition does not change the total momentum 
of the 67 system. 

Is such a change desirable? First note that after the 
P"]^ term is added, the photon momentum P^y^ is no 
longer parallel to its free-space momentum k. That is, 
the interacting photon is off shell. This unusual feature 
must be kept in mind if the modified operator is used. 

Secondly, the two types (nondynamical and dynam- 
ical) of electron momentum operators satisfy differ- 
ent commutation relations. Both operators satisfy the 
momentum-position commutation relations 



[p'^ -eA^{x),x,] 



(23) 



if A^(x) is momentum independent. However, between 
the momentum components themselves 



= ^||(VA||) + (V • E||)A|| - V^(£;i|A||) 



(17) 



[p^ - eAt'{x),p'' - eA^ix)] = ^eF'^^ 



(24) 
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the commutators vanish only for the field- free A'^^{x) 
that appears in the nondynamical momentum PJj^j. In 
contrast, not all three components of the dynamical mo- 
mentum Pjyjj are simultaneously observable jsj. These 
results remain unchanged in QCD. 

Thirdly, an experimental momentum measurement of 
the e or 7 bound in an atomic state requires the use of 
an external probe to eject the measured particle from the 
atom and put it on the energy shell for detection at a dis- 
tant detector. The detector cannot tell the past dynam- 
ical history of one detected free particle. It is only after 
extended measurements of suitable distribution functions 
when the past dynamical contents of the detected free 
particles can be deduced. In other words, unlike classical 
particle dynamics, the dynamics contained in a quantum 
wave function is not localized at one point in space but 
is spread out all over space in a probability density dis- 
tribution. 

Theoretical analysis requires calculations of relevant 
transition matrix elements where the full dynamics con- 
tained in the atomic states can in principle be included 
through their wave functions. In practice, at some stage 
one is often forced to use perturbation theory where the 
unperturbed states are often constructed from free-space 
operators. After that, the 67 interaction is included ex- 
plicitly order by order. In covariant perturbation the- 
ory, for example, internal as well as external lines are 
labeled by free-space 4-momenta. All ej interactions ap- 
pear as vertices where free-space energy-momentum con- 
servation is imposed. The propagators describing inter- 
nal lines also involve these free-space momenta without 
any interaction. Hence the nondynamical canonical mo- 
mentum operators are well matched to these theoretical 
and experimental circumstances. 

The interaction momentum term in the dynami- 
cal e and 7 momentum operators involve both particles. 
If we measure only one of the two particles, we do not 
know if the last vertex the particle comes from is not from 
other dynamical correlations in the atomic states rather 
than from this interaction momentum operator. If we 
measure both particles, we do not know if they actually 
come from the same interaction vertex, a knowledge that 
would qualify that vertex as the interaction momentum 
operator of interest. Moreover, perturbation theory gen- 
erally does not calculate both members of a canceling 
action-reaction pair of terms, but simply set their total 
contribution to zero, as is done explicitly in the nondy- 
namical canonical method. 

Thus it appears that our nondynamical momenta are 
simpler and perhaps more natural and more practical 
than the dynamical operators advocated by Ji ^ and 
Wakamatsu 



momentum operators of the last subsection. Invariance 
of the QED Lagrangian density C of Eq. under 
rotation gives rise to a certain conserved moment tensor 
TVf "^T whose spatial component M°^'^ integrates to the 
three-dimensional total AM 



J = S'^ + L'' + S'^ -t- 



E^{iLA^)d^x, (25) 



where = {i/2)ejmnl™'7" ^-nd L = xxp is the quantum 
orbital AM momentum operator. The photon operators 
can be expressed as (with details given in the Appendix) 



T 7 _L T '>' 
phys ' -^puroi 



S7, = / E_L X A_L d^x, 



T 7 

phys 



pure 



E^, (ilM ) d^x, 



V't(x X eA||)?/'rf^x. (26) 



We get two physical, gauge-invariant terms Sj| and Lpj^y^, 
and a nondynamical gauge-dependent term L^^.^, to be 
added to the electron orbital AM to make it manifestly 
gauge invariant: 



+ Lpure 



?/;^[x X (p- eA||)]^d^x. (27) 



The separation of Sjj' from Lp^yg has actually been known 
for many years (see, for example, [H, [H, |2^), as well as 
the implicit gauge independence of Coulomb-gauge quan- 
tization [11] • 

Ji's dynamical electron AM operator 

Ldyn = J {p- eA)ip]d^x 



_ T e r ej 
— ^nd - ^ ■ 



(28) 



is obtained by adding to our nondynamical li^^ a term 
—IJ''* that comes from the 67 interaction, where 



^67 



'^'''(x X eA±_)ip d^x. 



(29) 



The canceling partner L"^''' of this action-reaction pair of 
terms is then added to the photon momentum to give a 
dynamical photon AM operator 



C. Angular momentum operators in QED 

The construction of gauge invariant AM operators in 
QED proceeds in the same way as for gauge invariant 



dyn piiys 



(30) 



Equation (|26p shows that in momentum space, the spin 
operator Sj = 5|j'ek is built up entirely of 2-dimensional 
photon fields on the physical plane. This intrinsic photon 



7 



structure contains the hehcity operator Sj whose two 

eigenstates |e±) — tK^i ± 162) /\/2) lie on the physical 
plane. The helicity spin state jea) parallel to k is absent 
because the two transverse spin operators on the plane 
that can generate it are absent, i.e., = 0. These 
properties are caused by the fact that the massless photon 
has no rest frame. Only the realizable nonzero photon 
helicities are physically observable. 

On the other hand, the photon orbital AM L'' has all 
three intact components for any choice of the quantiza- 
tion axis . It is an external 3-dimensional operator act- 
ing on the photon fields whose directions are restricted to 
the physical plane. Normally a gauge invariant quantity 
like h'' can be expected to be physically observable. We 
shall show in the next subsection, however, that with all 
longitudinal vector spherical harmonics perpendicular to 
the physical plane absent, the orbital AM of a free pho- 
ton is only partially observable from the perspective of 
quantum measurement. 



D. Free photon operators and states 

The AM properties of a free photon become more 
transparent after canonical quantization in the Coulomb 
gauge d [H, m m when the vector potential 

takes the simple form 

A-l(-) = E (^•'^"^'■^ + ^^e''-^) ■ (31) 

k 

Here — (i,x), and the 3- momentum state labels k 
may take a non-Cartesian form to match the geometry 
of the quantization volume V. For the proper treatment 
of rotational properties, for example, it is best to use a 
spherical box. Then the AM eigenvalues contained in the 
k labels remain discrete as V ^ 00. (The same generic k 
notation can be used for a cubic or cylindrical box.) The 
transverse operators 

ak = E «kAeA (32) 

and destroy and create a single photon, respectively, in 
the Heisenberg (time-independent) representation. The 
time dependence appears solely in the four scalar product 
k ■ X ^ ujt — h ■ Ti, with uj — |k| for one-photon states. 

The fields E±{x) and B(a::) for free photons can be 
calculated directly from A±{x). Other operators can 
then be reduced easily to second quantized form. Thus 
the Hamiltonian, without the divergent zero-point energy 
term, and the momentum can be expressed in terms of 



number operators: 

H-< ^ ^ f (Ei+B2)d3a; = ^Lj7Vk, (33) 

k 

= J2^N^ = T. + A^k- ) , (34) 

k k 

iVk = a[-ak -A^k+ + Afk-. (35) 

The spin operator shows a characteristic spinning or 
vector product structure that can be simplified into two 
distinct final forms, a helicity form if ~ ek and a 3- 
dimensional form if is any direction: 

Sj[ - -iE4xak = Eek(A^k+-A^k-), (36) 

k k 

= EE(4)™^'""(^k)n = E^k^^k- (37) 
k m,n k 

Here m, n — 1,2,3 are Cartesian indices and (ak)n — 
Cartesian component in the chosen frame. 
The 3-dimensional form contains an expected quantum 
spin vector for spin 1 particles, 

S = ^ejSj, where {Sj)mn = -if^jmn, (38) 

j 

that is a 3 X 3 matrix in spin space with Cartesian compo- 
nents Sj that are 3x3 Hcrmitian matrices defined by the 
Levi-Civita symbols from the cross product. It satisfies 
the standard quantum AM algebra S x S = iS. In this 
matrix notation, the 3-dimensional vectors ak,aj^ in the 
last expression of Eq. ((37|) are column and row vectors, 
respectively, in spin space. 

Finally, the orbital AM operator has the same general 
structure: 

= E E '^kAL flkA = E 4Lak. (39) 

k A=± k 

Note the subtle difference between the expression for L,'* 
where L = — k x x acts on the components akA of the 
2-vector ak on the physical plane, and S''' where S acts 
on the unit vectors contained in the same vector ak- 
The last term in Eq. can also be considered a matrix 
expression with L a non-matrix quantity, or a diagonal 
matrix in spin space if an optional 3x3 identity matrix 
is added. The 3-dimensional forms for both 'L'^ and S'' 
are needed in real space where the physical plane cannot 
in general be isolated explicitly. 

Three of the operators, H'*,P'* and Sj|, are functions 
only of the number operators that commute among them- 
selves. Hence the n-photon number state is also 
the simultaneous eigenstate of these operators with the 
respective eigenvalues 

n 

S-E^'' K = ^k„ A = E^^e,, (40) 

i=l i i 
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where = of the i-th free photon. In particular, the 
one-photon state of momentum k can be written in the 
hehcity form (23 | 



= k\)=al\0). 



(41) 



This photon state carries only 4 labels or constants of 
motion: , j — 1,2,3 (including uj = |k|), and a single 
helicity. The Cartesian components of the spin vector Sj 
in momentum space commute among themselves j26l |. 



S7, ,5; ] =0, 

1 1 m ' 1 1 n -I ' 



(42) 



where Su 

\\r 



Sj itself is a generator of 



two-dimensional, Abelian rotations about et- It is not 
a normal 3-dimensional spin operator in momentum or 
real space. 

Helicity states of good total angular momentum for a 
free photon can also be constructed |24] : 



\kJMX) = 



2J + 1 
in 



where k = ek, and -Ojj'J^ is a representation matrix ele- 
ment of the finite rotation operator. Note that there are 
also four state labels {k, J, M, A) here, and that a state 
of good J does not have a preferred momentum direction 
Gk or k. 

All these well known expressions take such simple 
forms only because the photon is free. If the 67 interac- 
tion is also included, these expressions will become much 
more complicated. 



d2ki5i^i(k)*|kA), 



(43) 



E. Photon multipole radiation 



The helicity plane- wave state \a 
is also an eigenstate of L7. = Gk 



T 7 

phys 



kA) of a free photon 
with eigenvalue 
Ml = 0, because the state is axially symmetric about 
k. As such, it is one member of a set of states in cir- 
cular coordinates with any integer eigenvalue Ml called 
Laguerre-Gauss modes in optics ^29,,^]. The modes with 
Ml ^ are physically observable helical waves [3l| that 
do not concern us in this paper, but they do reinforce our 
understanding that the component Lj of the orbital AM 
of a free photon along Gk is measurable. 

What about the remaining components of L'''? The is- 
sue can be addressed readily by using the elegant method 
of Berestetskii [s^l for constructing coupled AM states 
\{LS)JMj) or vector spherical harmonics that include 
the multipole radiation of massless free photons. The 
TO-th component of any vector field, say in momentum 
space, can be written in the form 



(44) 



where a vector component Vm can be separated from the 
remaining scalar field /(k) common to all components of 



A.. Under rotation, the vector V behaves just like the 
momentum k. In fact, it can be taken to be k itself. The 
vector k does not commute with the quantum orbital AM 
operator L = — k x x, but satisfies the commutation rela- 
tions [Lj, k„i\ — isjmnkn- Similar commutators therefore 
must hold for any vector V: 



(45) 



(The commutator is also valid in real space where the 
position vector plays the role of k for the vector field. 
We use the same font for a quantity in both momentum 
and real spaces when its properties are very similar in 
both spaces.) 

The photon orbital AM density appearing in the third 
equation of Eq. (PBl) then has the structure 

£*J.{V^f) = C(KnL + Gjiej™„K)/. (46) 

In matrix form, it reads 

£^t{Yf) = 5t(VL-SV)/, (47) 

where V is a column vector and 5^ is a row vector. One 
can then extract from it the Berestetskii identity as the 
operator equation in momentum or real space 



J(V/) = VL/, 



(48) 



where J = L -f S is the quantum total AM operator. All 
these quantum AM operators satisfy the standard AM 
algebra: L x L = iL, S x S = zS, and J x J = zJ. 

For the special case where the scalar field / is the 
spherical harmonic i^/M, one finds 

3^{YYjm) = VVYjm = J{J+1){YYjm), 
iz{YYjM) = YUYjm = M{\Yjm). (49) 

Hence the vector state |V1^7m) is a coupled \{LS)JM) 
state. The method actually works for any vector particle, 
massive as well as massless. For massive particles, one 
can show that the result is identical to that obtained from 
Clebsch-Gordan coefficients '34] . 

For massless photons, however, the vector fields must 
reside on the physical plane. Hence the out-of-plane lon- 
gitudinal vector fields kYjM are nonphysical and unreal- 
izable. The permissible fields are two types of transverse 
fields often taken to be LIjm and Gx^jm. A transverse 
field of the first type is called a magnetic or MJ radi- 
ation, and may be denoted by the spectroscopic symbol 
'^{L = J) J, with a single L = J value for each J. A trans- 
verse field of the second type is an EJ radiation made up 
of a specific linear combination of the two states ^{J—l)j 
and ^{J + l)j of different L values, a combination that 
is orthogonal to the linear combination found in the lon- 
gitudinal field of the same J. 

The unavailability of the longitudinal field for free pho- 
ton states means that it is not possible to construct the 
uncoupled \LMl)\SMs) states from just the two trans- 
verse physical photon states. That is why one cannot 



9 



measure simultaneously and exactly both L and Ml val- 
ues of the photon orbital AM in free space. The photon 
Ml is known to be zero along the momentum direction 
k. The value L = J for MJ radiations is known if J is 
known. On the other hand, each EJ radiation contains 
two possible L values. The average values of and Ml 
can be calculated if J, M are known, but not their exact 
values. So the orbital AM of a photon emitted by an 
atom into free space is only partially observable. 

More generally, there can be other complications to 
the quantum observability of the orbital AM even for 
massive matter particles such as electrons and nucleons 
when they are inside composite systems. For example, if 
L does not commute with the Hamiltonian, as happens 
in the deuteron when a tensor nucleon-nucleon force is 
present, the relative orbital AM L of an energy eigen- 
state such as the deuteron ground state is not a constant 
of motion (characterized by simultaneous quantum num- 
bers) of the quantum energy eigenstate. For this reason, 
the deuteron D-state probability po = (L^)/6?i^ is not 
the eigenvalue of a physical observable that commutes 
with the Hamiltonian. However, it can be indirectly de- 
duced if somehow we can measure (L^). 

A final remark: For quantization in a covariant gauge, 
the longitudinal and time-like photons are both unphysi- 
cal and not directly observable as on-shell photons. How- 
ever, they do give rise to observable effects such as 
the Coulomb interaction among static charges [sl] and 
Coulomb excitations/deexcitations of atomic and nuclear 
states. However, these effects are also correctly described 
for quantization in the Coulomb gauge using only trans- 
verse photons. 



III. GAUGE INVARIANT OPERATORS IN QCD 

In QED, the gauge dependent and nondynamical 
part ^iid of the vector potential A'^ resides in the 2- 
dimensional energy-momentum (03) subspace orthogo- 
nal to the physical (12)-plane that contains the gauge- 
independent part A.±. In QCD, one intuitively expects 
the same separation because color is an internal attribute 
independent of spacetime itself. We obtain in Sect. A an 
explicit perturbative solution for the nondynamical And 
and find that it is indeed confined to the (03)-subspace. 
In Sect. B, both nondynamical and dynamical momen- 
tum operators are found to have appearances similar to 
the QED operators. However, the color electric field in 
the gluon operators differs from the QED field in an ex- 
tra term quadratic in A and proportional to the QCD 
coupling constant g. This term describes an interaction 
that changes the value of the gluon momentum itself. 
The same nonlinear electric field term is found in Sect. 
C to add nonlinear interaction structures to our orbital 
AM and total AM operators in QCD that are otherwise 
nondynamical. 



A. Nondynamical vector potentials in classical CD 

The color Maxwell field is F^"' = F^^Ta, where 

Fr = O'^A: - a^A^: + gCabcA^^A-^^. (50) 

The color operators Ta are 3x3 matrices in color space 
that satisfy the SU3 commutation relations 



(51) 



where the structure constant Cabc is totally antisymmet- 
ric in its indices. The nondynamical (nd) vector potential 
= ^nda-^a gluou field is defined by the condi- 

tion that it generates no color field: 

KL = d^K - + gCabcA^A: = 0, (52) 

where the subscript "nd" has been dropped from for 
notational simplicity. We do not know any analytic so- 
lution of this nonlinear differential equation for j4J^, but 
look instead for a perturbative solution in powers of the 
QCD coupling constant g: 



00 



n=0 



(53) 



The first term A^a'^^ is just the linear QED solution, be- 
cause it satisfies the QED zero-field condition d^A'^ — 
d'^A'^ = 0. The remaining terms are solutions of the 
differential equations 

n~l 

Qt^Ain)^ _ g^^(^n)M ^ (J^^^ 4''^^4"-l-J>. (54) 

In the Fourier space A:^ = (w, 0, 0, |k|), Ea. ([5i|) reduces 
to the matrix equation 



for 71 > 0. Here 

n-l 



(55) 



X 4"-i-J>(fc')rf*fc'. (56) 

have been obtained from the convolution theorem. Note 
that £1°)^" = Ti:^^ = 0. 



Equation (|55p equates two 4x4 antisymmetric matri- 
ces. On the right, the matrix 7?.i" will be shown 
inductively to have only two nonzero matrix elements, at 
positions fiv — 30, 03 where 63 = e^. On the left, the 
matrix has the structure of a general antisymmetric 
field tensor with 6 zero matrix elements, 4 on the diagonal 
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and 2 at positions 12 and 21 (because = k'^ = 0). The 
remaining 10 nontrivial matrix elements are in 5 anti- 
symmetric pairs, 3 electric field components and 2 mag- 
netic induction components. The component -B|| = B3 
vanishes in this QED-like field tensor because magnetic 
monopoles are absent in QED. 

With the matrices written in the Cartesian order 
(3, 0; 1, 2), the 4 matrix elements in one of the two 2x2 
off-diagonal block submatrices is assumed temporarily to 
satisfy the algebraic equations 



(57) 



where /3 



1,2, (7 ~ 3,0. This assumption, that 
0, will be justified by induction later. The 
unique solutions of Eq. ([57)1 are then Aa^''^ = 0, for 
(3—1,2. In other words, the temporary assumption is 
that the driving term 7?.i" '^''^'^ does not have any block 
off-diagonal matrix element to drive aIT"^^ onto the phys- 
ical (12)-plane. In this way, the nondynamical Aa stays 
in the (03)-subspace orthogonal to the physical plane. 
The same solution can be obtained from the other off- 
diagonal block submatrix. 

The 5th and last equation is in the (30) subspace alone 
and reads 



(58) 



A second equation is needed to solve uniquely for the 
two unknowns in the n-th order. This is provided by 
choosing a gauge, but so far we have considered only 
the velocity gauge condition (|4]) for QED that is good in 
QCD only for order n = 0. For n > 0, one can in principle 
choose any velocity gauge in every order n in a totally 
random manner, but such a choice gives infinitely many 
possibilities that cannot even be written down explicitly. 
The simplest choice is the opposite extreme, namely, to 
use the same velocity gauge in every order of g: 



A(n)3 _ „[f^A{n)0 



(59) 



where s = ±1 and the gauge parameter a 
defined by the gauge velocity Vg. 

The complete zeroth-order nondynamical vector po- 
tential [in the component order (3, 0; 1, 2)] is known from 
QED to be: 



Ai°) = (^r,^(°)";0,0) 



(60) 



for every color a. The antisymmetric matrix TZli^ thus 
have two rows and two columns of zeros where the matrix 
index is on the physical plane, i.e., /i = 1, and 2. Hence 
it has only two nonzero matrix elements, at positions 30 
and 03. This verifies the assumed property of 7?.i"^ for 
n = 0. 

The first-order equation ((55)) (with n = 1) and the 
order-independent gauge condition (1591) can now be 



solved for the two unknowns. The result, 

Ikl 



(1)0 _ 



-4, 



|kp — sau. 



.77(0)30 



(61) 



is a vector confined to the (30)-plane, just like A!f K The 

antisymmetric matrix TZ'^^ constructed from it also van- 
ishes everywhere except at positions 30 and 03. The 
second-order solution can then be obtained in exactly 
the same way. We can thus bootstrap our way up the 
perturbative expansion to give the complete solution 



A". - 



A^. 



|kp — saui 

aO 

Ikl 



2 '^a 



30 



(62) 



n=0 



confined to the (30)-plane. The nonlinear dynamics 
changes the actual functional form on this plane, but 
preserves unchanged its avoidance of the physical (12)- 
plane. 

Any color vector potential in a chosen Lorentz frame 
can thus be written as A^^ — ^quan + ^nd' where ^JJuan is 
its dynamical part in the gauge where second quantiza- 
tion is realized. Here the color label is unimportant and 
is omitted. For canonical quantization in the Coulomb 
or physical gauge, we have A^^^^ = A'^^^^^ = A'^. With 
the same velocity gauge condition as QED, it follows that 



^phys 



Aj^ holds in any color. 



The use of the Coulomb gauge also simplifies the non- 
linear (NL) part of the field tensor. With A\\ = A^ = 
-^physNL is a 4 X 4 matrix that is not only antisymmet- 
ric, but is also zero in the row and column involving a 
j = 3 component. Hence it contains only 6 nonzero ma- 
trix elements or 3 nonzero field components: , and 
. The presence oi ^ Q means that the nonlin- 
ear color dynamics generates color magnetic monopoles. 
As a result, the CCD momentum density contains addi- 
tional terms including one Enl x Bnl = E^nl x B||nl 
that gives a contribution in a direction perpendicular to 
k. Similar complications arise in the AM operators. 

A final note: Since the differential equation ([5^ has 
a driving term quadratic in A^, it seems likely that the 
equation has another solution. However, this cannot be a 
perturbative solution in powers of g because the pertur- 
bative solution obtained here is the unique perturbative 
solution for the chosen gauge. If the second solution ex- 
ists, it will be a nonperturbative solution. 

While it is interesting that QED-like gauge conditions 
such as Eq. ([55]) also works in QCD, it is worth not- 
ing that the formalism developed here works for more 
complicated gauge relations as well. In the rest of the 
paper, we shall consider only the QED-like perturbative 
solutions obtained in this subsection. 
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B. Momentum operators in QCD 



The color Maxwell field ((50)) is conveniently written as 
two 3- vector fields 



(63) 



V X A„ - -gCabc (Ab X Ac) . (64) 



Of the color Maxwell equations ^^F^" + gCabcA^fiF^" 
, we shall only need the color Gauss law: 



V-E„ 



Pa + gCabcAb ■ Ec 



(65) 



where in QCD, pa = g4'''Taip. 

The derivation of momentum and AM operators in 
QCD differs from that already presented for QED only 
by the addition of nonlinear terms proportional to g. The 
similarity allows the separation = -f*phys + Ppuro to 
be made in any velocity gauge when quantization is re- 
alized on the physical (12)-plane perpendicular to et- 
Again P^u^c can be simplified by using the QCD ver- 
sion of the vector identity (|T7)) by considering the special 
case Bpuica = 0, for which Apmca = A||a and therefore 
both terms on the right of Eq. ([M]) vanish separately. 
Note that the nonlinear term on the right of Eq. is 
nonzero only when the spatial direction of the vector po- 
tential Af, depends on the color label b. We do not know 
if such solutions with space-color correlations exist. 

Using E — E||(j, we find with the implied sum over the 
color label a 

= £;||a(VA||J + (V • E|| JA||, - V-'"(i?|^A||,). (66) 

The color Gauss law (|55|) has an additional nonlinear 
term on the right, but its contribution to Eq. (1661) on 
the right is gCabcA\\ciA\\b^\\c — 0. Hence all the QED 
decompositions and rearrangements can be taken over 
wholesale by substitutions of the type eA^ gA^ = 
gA^gTa and by summing over an additional color label 
where appropriate. 

For example, the manifestly gauge-invariant, nondy- 
namical quark momentum operator takes the same form 
as Eq. (Ull) 



p9 ^ p9 

I pure 



ip\-p - gAMd^x 



(67) 



when the color content is not made explicit. However, the 
color structure has to be made explicit to display explicit 
gauge invariance: The quark field tp is & column vector 
in color space. Under a nondynamical gauge transfor- 
mation, it is changed to V' = Uijj, where U = exp (zA) 
is a 3 X 3 matrix in color space because A = XaTa- The 
quark momentum density "0^ {p^ — eA^^^^)%l} then remains 
unchanged if (p^ - eA^^ro)' = Uip'^ - eAf^^^^)!!"^ . Just as 
in QED, both fc'^ and ^pm-o have only nonzero /i = 0, 3 
components in momentum space. Their identical math- 
ematical structures explain why they appear naturally 



together. The density tp'^A±tp in the dynamical opera- 
tors of Ji 01 is gauge invariant if A^ = UA± [/^ . 

With P^uro removed, the remaining gluon momentum 
density £'^(VA^) — i?ij^(VA||^) can be written as three 
terms of spatial structures _L||,||_L, and _L_L. The two 
cross terms vanish separately, leaving 



phys 



ps 

physL 



pg 

physNL' 



p9 



C=L,NL- 



(68) 



The linear (L) and nonlinear (NL) terms come from the 
first and second terms, respectively, of the color electric 
field on the physical plane 



E±a — —dtA±a — gCabcAlA±c, 
= Ej_aL + E^aNL- 



(69) 



The NL term is an interaction term proportional to g 
that is not present in QED. It contains the nonlinear 
momentum density in momentum space of 



gCafccK * -4lJ(k)[zk^iJk)] ^ 0. 



(70) 



If second quantization is still defined by A of Eq. (PT|) 
for each color, E_lq remains linear in a and a^^, but the 
interaction contained in the nonlinear term gives rise to 
gluons of arbitrary momentum k' inside the convolution. 
[The convolution contributes to Eq. (fTO)) except when 
k' = k.] Hence gluon momenta are more complicated 
than photon momenta. 



C. Angular momentum operators in QCD 

The construction of AM operators proceeds similarly 
to QED. The details will also be given in the Appendix. 
The results are summarized here. 

The expression obtained directly from the QCD La- 
grangian is 



E, X A„ d^x + j Ei{{LA^^) d^x, (71) 
The gluon operators can be re-written as 

^ — + %hys + ■'^purei 

Sff = / Ei, X A^a d\, 
L9_„ = - / 7/'t(x X gA||)V'd^x. (72) 



T 9 _ 

phys 



pure 



We get two physical, gauge-invariant terms Sy and Lpj^^^,, 
and a nondynamical gauge-dependent term to be added 
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to the quark orbital AM to make it manifestly gauge 
invariant; 

^nd - + -Lpure 

= y^/^tx X (p-.gA„)i/;d3x. (73) 

Ji's dynamical quark and gluon AM operators can be 
obtained from the electron and photon operators by a 
simple change of names when the color index is not shown 
explicitly. 

Returning to our nondynamical gluon operators (leav- 
ing out the pure gauge term L^j^^g for simplicity because 
it is not involved in the following discussion), we note 
that the following results for their nonlinear terms: 

S||NL = -gCabcJ (AlA^,) X A^a X = 0, (74) 

LphysNL = ~9Ca,cJ{A'^AiJ^LAi^d^x^0. (75) 

That is, S^j^L Eq. ([71]) vanishes for all solutions A±a 
whose direction in real space x is the same for all colors. 
In Eq. (|75p . the color sum would have given nothing if 
it were of the form CabcXbXc — with two identical and 
equivalent vectors Xa and Xh that can be interchanged 
freely by re-labeling the color indices. However, the ac- 
tual expression involves A-'j_^ and LA^^ that are not iden- 
tical vectors. It is therefore nonzero. 

A more general result can be derived for the sum 
JphysNL = LphysNL + S^nl of Eqs. d ESJ. The vector 
form of the Berestetskii identity Eq. (|A4[) can be written 
in real space for A_La(x) = V_L/a(x), where /o(x) is a 
scalar function: 

Eia'UVl fa) + Ei„ X A^a = {E^a ' V^)ltfa. (76) 

Hence 

JphysNL = -gCabcJiA°Utfa)^^ld'x^O, (77) 

because fc and L/^ are different functions in general. For 
example, if / = ^lm, then L/ can be written as a sum 
of three terms, each proportional to e^Y^Af, e+Y^AZ-i 
and e_Y£Af_|_i, respectively. The first term contributes 
nothing to Eq. (|77p , but the remaining two terms do con- 
tribute. Thus Jphys ^1^*^ has a more complicated struc- 
ture than Jpiiys in QED. 

IV. CONCLUDING REMARKS 

Ji in his published Comment [TtI on our recent pa- 
per |9j on angular momentum (AM) operators in gauge 
theories has expressed a number of concerns. First, he 
claims that our "new notion of gauge invariance clashes 
with locality and Lorentz symmetry." Let us first re- 
spond to his concern about Lorentz symmetry. As we 



have explained in the Introduction, manifest Lorentz co- 
variance and gauge invariance can coexist explicitly only 
in a covariant gauge such as the Lorenz gauge. For quan- 
tization in a non-covariant gauge, we need to retrace our 
steps by transforming to a covariant gauge in the chosen 
Lorentz frame and use the latter's Lorentz covariance to 
restore explicitly Lorentz symmetry. The point is that 
once Lorentz covariance is established initially, it does 
not have to be maintained in subsequent steps. So there 
is no need for us to look at the vector potential in an- 
other Lorentz frame when we quantize a gauge-invariant 
theory in a non-covariant gauge. 

However, we find hidden Lorentz symmetries in our for- 
malism in Sect. Ill Al First, is a Lorentz scalar, where 
A± is the gauge-independent part of any vector poten- 
tial A in our chosen family of gauges in 4-momentum 
space, the part that lies on the 2-dimensional physical 
plane perpendicular to the photon/gluon momentum k. 
This is an interesting and desirable feature, even though 
it does not involve the whole vector potential. 

Ji's concerns about both Lorentz symmetry and local- 
ity are also resolved by our choice to second quantize 
only on the physical plane. We have shown explicitly in 
Sect. Ill DI that for any choice of gauge among our chosen 
family, second quantization nevertheless proceeds like a 
standard canonical quantization in the Coulomb gauge 
for the gauge- independent part ^phys = •A±- No un- 
usual or unacceptable nonlocal operators appear. The 
simple Coulomb-gauge structure is always obtained be- 
cause we do not use the complete vector potential A'^ 
directly. Ji's concern may apply to quantization for the 
complete A^, a procedure we do not recommend. More- 
over, our quantization procedure is mathematically iden- 
tical to Coulomb-gauge quantization and shares the lat- 
ter's implicit Lorentz-invariant character in the resulting 
quantum field theories 22, 23]. 

Thirdly, we have taken considerable space in the Intro- 
duction to explain why the nondynamical gauge transfor- 
mation in the vector potential and the associated path- 
independent phase factor in the fermion field contain the 
essence of the concept of a gauge invariance that does not 
change the assumed dynamics of the problem. After this 
demonstration of gauge invariance, calculations can then 
be made in any convenient gauge. In Sect. Ill Bl we have 
pointed out some theoretical and practical advantages of 
not including dynamics in the definition of operators of 
certain wholly or partially observable physical attributes, 
at least in QED. 

However, there is a different gauge principle, one that 
gives rise to the local gauge theory of interactions, where 
dynamics can be generated by the use of dynamical, path- 
dependent or nonintegrable phase factors that multiply 
into a dynamics-free fermion expression. From our per- 
spective, nothing forbids the use of the dynamical oper- 
ators advocated by Ji and others. We have described in 
this paper, especially Sect. HE, the advantages of our 
nondynamical operators based on p^^ — gAI^^^^ over Ji's 
dynamical operators based on p^^ — gA^ . Even the math- 
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ematical structure of our operators is simpler because 
in 4-monientuni space, both fc^ and ^pm-e have nonzero 
components only for ^ = 0, 3 (or ||), i.e., in the subspace 
orthogonal to the physical plane. Which type of oper- 
ators will be more convenient to use will be settled by 
future practices. 

We have also pointed out in Sect. Ill Bl an essential 
difference between classical and quantum mechanics in 
their dynamical descriptions. The dynamical momentum 
appears naturally in classical particle dynamics where a 
point particle carries around its full dynamics as it moves 
in space. In contrast, an observed quantum particle is 
nondynamical in the sense that it is on the energy shell 
and does not contain any information on its past dy- 
namical history in the quantum source system where it 
existed before the detection. One has to measure many 
such particles from the same quantum source in order 
to deduce their common dynamics from their measured 
distribution functions. 

Among the issues under popular discussion is the ques- 
tion of whether the orbital AM and spin of photons and 
gluons are separately observable. For direct measure- 
ments on free photons, the answer is unambiguous, as 
we have described in Sect. Ill El The transverse pho- 
ton states can be expressed in terms of one or two states 
of the type \{LS)JM), where S", J, M are good quantum 
numbers. In MJ photons, only one state L = J appears, 
whereas in EJ photons, two states with L — J ±1 appear 
in a known linear combination. Hence S, J, M are com- 
pletely observable, and L is observable for MJ photons 
and only partially observable for EJ photons. In these 
\{LS)JM) states. Ml or Ms is not separately observ- 
able, but the combination M = AIl + Ms is observable. 
However, in states of good momentum k, M^ — and 
Ms = A = ±1 along et are separately observable. 

For gluons, however, our analysis shows that even 
our nominally nondynamical P and J operators are not 
free of interaction contributions. This is because the 
color electric field E_l appearing in them contains a term 
quadratic in the vector potential A'' and proportional to 
the color coupling constant g, as shown explicitly in Eq. 
([M)) . The fact that gluons are always confined and never 
free can cause additional complications. The present sta- 
tus of, and open problems in, the experimental measure- 
ment of quark and gluon contributions to the proton spin 
have been reviewed recently by Bass [35j . 

Finally, Ji [l3| has pointed out the energy scale de- 
pendence of various quantities involved in the nucleon 
spin problem. This dependence comes from the fact 
that all physics, including the strong coupling constant 
q;s(/x) = (7^(/i)/47r, is a function of the energy scale ^ at 
which the measurement is made. He has emphasized that 
even the quark spin, which has gluon contributions from 
the Adler-Bell- Jackiw triangle anomaly, is not free of this 
complication. Indeed, as discussed in the last paragraph, 
all QCD momentum and angular momentum operators 
are (7-dependent. The lesson here seems to be that the- 
oretical analyses should match the energy scale of the 



experimental data. However, various commutation rela- 
tions that simplify calculations at one energy scale can 
still be used in spite of their intrinsic scale dependence. 
Again, the question will be settled by future practices. 

In contrast, the gauge invariance addressed in this pa- 
per is a relatively simple issue. The simple answer given 
here is that this gauge invariance is actually a nondy- 
namical concept in both QED and QCD. 



Appendix A: Angular momentum operators in QED 
and QCD 

In this Appendix, we show how the AM operators cal- 
culated from the Lagrangian in QED and QCD are re- 
duced to a form suitable for canonical quantization in the 
Coulomb gauge. 

For QED, Eq. ([25)1 is obtained from its parent expres- 
sion (|25p if their difference (expressible as two groups of 
terms) vanishes. The first group of terms is 

+ j Ex A±d^x = 0. (Al) 

This equation can be derived by using the differential 
vector identity pT|) for the special case B = Bpurc = 
B|| = 0, A = A||: 

= E^{x X V)A| + (V • E)(x X A||) 

+ E X A|| - V^' [E^ (x X A|| )] . (A2) 

The spatial integral of the last term on the right vanishes 
when it is changed into a surface integral at 00 where the 
integrand vanishes. Equation (jAll) then follows with the 
help of the Gauss law 'V ■ E — p — eip'^ip. 

The second group of terms is more transparent when 
written as a momentum integral 

J [SfiiLA^^) + £lx Ai_] d^k = 0. (A3) 

To verify its vanishing value, we write the Berestetskii 
identity (H71) for two arbitrary vectors £ and A = V/, 
where / — /(k) is a scalar field, in the alternative vector 
form 

Ei*t{V^f) = {£* ■ V)L/ + i£* X A. (A4) 

For the special case £ = £\\, A — Ajl, the identity shows 
that the integrand in Eq. (jA3p is equal to the 5| • Vj_ 
term which vanishes at every value of k. This completes 
the derivation of Eq. (gl]) from Eq. (^5]) . 

For QCD, the second group of terms [the left side of Eq. 
(|A3p ] also contribute nothing because Eq. (jA4p holds for 
each color. For the first group of terms [the left side of 
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Eq. (IA1[) ]. we need to use the differential vector identity 
(|A2[) with an imphed sum over the color index a: 

= i?^(xx VM|^ + (V.E,)(xx A||J 

+ E, xA[|,-VJ"[i?^.(xx A||J]. (A5) 

The only real difference from the linear QED structure 
is the additional contribution coming from the nonlin- 
ear term in the color Gauss law (|55|). With the implicit 



sum over the color index a, this nonlinear contribution 
vanishes because 

gCafcc(A||6 • Ec)(x X A||a) 

= ge^ej^sx^'iCabcAlAlE^J = 0, (A6) 

where e|| — e^. Hence the first group of terms vanishes 
also for QCD. Thus Eq. ^ follows from Eq. dZJ)- 
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